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Smeary 

The  mean  residual  life  (mrl)  function  tp(z)  gives  the  expected  remaining  life  of  a 
patient  at  age  x,  where  the  underlying  failure  distribution  is  F.  ReliabiUste  and 
bicmetrioiane  have  found  it  ueefitl  to  categories  failure  distributions  by  monotonioity 
properties  of  the  mrl.  In  particular  F  is  said  to  be  a  decreasing  mean  residual  life 
(dmel)  distribution  if  t^fO)  is  finite  and  for  all  0  4  e  4  t,  tp(e)  i  tp( t) .  If  the 
preceding  inequality  is  reversed,  F  is  said  to  be  an  increasing  mean  residual  life 
(imrl)  distribution .  Hollander  and  Prosohan  (1875)  have  derived  tests  of 
(1)  Hq:  F  is  exponential,  versus  H^:  F  is  dmrl,  and  (2)  HQ  versus  B F  is  imrl.  Their 
tests  are  based  on  a  (complete)  random  eomple  Xj,  ...,  Xn  from  F.  Often,  however,  data 
are  incomplete  due  to  withdrawals  from  the  study  and  patients  still  surviving  at  the  time 
the  data  are  analysed.  In  this  paper,  we  generalise  the  Hollander-Proeohan  tests  to 
accomodate  randomly  censored  data.  Asymptotically  distribution-free  teste  of  Hq  versus 
Hj  and  Hq  versus  H£  are  provided.  The  efficiency  lose,  due  to  the  presence  of  censoring, 
is  also  investigated. 


1.  Introduction 

Let  X  be  a  randon  variable  denoting  the  tine  to  the  occurrence  of  an  end-point  event. 
The  distribution  function  (df)  F(x)  ■  P(X  s  x)  is  called  the  failure  distribution  of  X; 
Rp(x)  ■  1  -  F(x)  is  the  survival  function.  The  nean  residual  life  («rl)  function 
corresponding  to  F(x)  is 


*p(x)  ■  { £Rp(u)du}/Rp(x) . 


(1) 
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t  ♦ 


The  quantity  cp(x)  is  the  ejected  remaining  life  of  the  patient  at  age  x,  and  thus 
ccxild  also  be  written  as  ep(x)  *  E(X  -  x|X  >  x) . 

The  mrl  is  a  function  of  central  importance  in  reliability,  biometry,  and  other 
disciplines  where  survivorship  studies  arise.  Some  examples  where  knowledge  concerning 
Cp(x)  is  desired  include  the  following. 

(i)  An  engineer  is  testing  a  component  to  study  its  failure  distribution.  Natural 
questions  include:  what  is  the  expected  remaining  life  of  the  coeponent  given  that  it 
has  functioned  properly  for  10  days?;  for  100  days?;  and  so  forth. 

(ii)  A  biometrician  is  studying  the  survival  data  of  patients  receiving  estrogen 
for  treatment  of  prostate  cancer.  The  biometrician,  physicians,  and  patients  are  interested 
in  the  mean  residual  life  at  12  months,  at  18  months,  at  24  months,  etc. 

(iii)  A  sociological  experiment  is  performed  in  which  the  variable  under  study  is  the 
tine  X  until  a  victim  is  extricated  from  a  'taishap"  by  receiving  help  from  a  passerby.  How 
long  will  a  victim,  on  the  average,  have  to  wait  if  he  has  received  no  help  in  the  first 
10  minutes?,  twenty  minutes?,  etc. 

Obviously  many  other  examples  where  study  of  the  mrl  arises  can  be  cited,  and  the 
above  examples  can  be  rephrased  in  different  contexts. 

Estimation  of  the  mrl  has  received  considerable  attention  in  the  literature. 

Sippose  a  random  sample  Xj,  ...,  X|1  is  available  from  F.  A  natural  nonparametric 
estimate  of  ep(x)  is 

m 

Spto  “  CCJ^OOduJ/RpOOJJ  •  ICx  <  X(n)J,  (2) 

where  X^  •  maximum  CXj,  X2,  ...»  X^J,  ICx  <  X^j  ■  1  if  x  <  X^  and  0  otherwise,  and 
RpCx)  is  the  empirical  survival  function  of  the  X's.  That  is,  Sp(x)  -  n"1.  (mater  of  X’s 
in  the  sople  >  x). 


Bryson  and  Siddiqui  (1969)  plot  ep  for  survival  data  of  leidtamia-petients.  Yang 
(1978)  established  strong  consistency  of  %  on  a  finite  interval  CO,  TJ  and  also  showed 

L  w 

that  the  associated  mrl  process  n  (ep(x)  -  tp(x)}  converges  weakly  to  a  Guassian  process. 
Hall  and  Nellner  (1979)  strengthen  Yang's  results;  in  particular  they  extend  her  weak 
convergence  result  to  the  positive  real  line.  Hall  and  Wellner  also  4erive-nonpara*etric 
sinultaneous  confidence  bands  for  ep(x).  Chiang  (1960)  focused  on  estimating  the 
discretized  life- table  version  of  the  mrl  function. 

Testing  for  properties  of  the  mrl  has  received  less  attention  in  the  literature  even 
though  reliabilists  have  found  it  useflfr  to  categorize  failure  distributions  according  to 
monotonicity  properties  of  the  mrl . 

A  failure  distribution  F  is  said  to  be  a  dsereasing  mean  residual  life  (<karl) 
distribution  if  the  mean.Cp(O),  is  finite  and, 

Cp(s)  a  Cp(t) t  for  all  0  s  s  s  t.  (2 

F  is  said  to  be  an  increasing  mean  residual  life  (imrl)  distribution  if  ep(0)  is  finite 
and  the  inequality  in  (3)  is  reversed. 

The  class  of  dmrl  distributions  is  useful  for  modeling  situations  where  lifelengths 
c<!  items  deteriorate  with  age,  whereas  the  class  of  imrl  distributions  is  appropriate  for 
models  where  lifelengths  improve  with  age.  The  boundary  mashers  of  the  dmrl  and  imrl 
classes,  obtained  by  requiring  equality  in  (3),  are  the  exponential  distributions.  Of 
course,  the  exponential  distributions  can  be  used  to  model  situations  where  lifelengths 
neither  improve  nor  deteriorate  with  age.  Equivalently,  they  are  the  distributions  for 
which  the  mrl  function  ep(x)  is  constant  for  all  x. 

Hollander  and  Proschan  (1975)  C»»(1975)]  consider  the  problem  of  testing 

Hp:  F(x)  »1  -  axp(-x/u),  x  a  0,  a  >  0,  w  unspecified,  (4) 

versus 


H^:  F  is  a  dmrl  distribution  and  is  not  exponential, 


(5) 
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using  «  ranks  sample  Xj,  ...,3^  from  F.  Significantly  large  values  of  their  test 
statistic  V*  lead  to  the  rejection  of  Hq  in  favor  of  H^.  Significantly  small  values  of 
V*  lead  to  the  rejection  of  Hq  in  favor  of 

H{:  F  is  an  iarl  distribution  and  is  not  exponential.  (6) 

HP(1975)  establish  asymptotic  normality  of  V*  and  also  provide  critical  values ,  obtained 
by  Monte  Carlo  simulation,  for  n  -  2(1)20(5)50.  Exact  critical  values  for  n  -  2(1)20(5)60 
are  given  by  Langenberg  and  Srinivasan  (1979) . 

Bryson  (1974)  has  suggested  a  test  of  Hq  vs.  H£.  However  he  gives  critical  values  only 
f or  n  -  10,  15,  20,  25,  and  30,  he  does  not  derive  the  asymptotic  distribution  of  his  test 
statistic,  and  his  test  can  not  be  easily  converted  to  a  test  of  Hq  versus  . 

In  this  paper  we  generalize  the  HP (1975)  test  to  accommodate  the  randomly  censored 

model.  Thus,  instead  of  a  complete  sample  3Cj . X^,  we  are  only  able  to  observe  the 

pairs 

Z.  -  minOtp  Tj),  -  1  if  Z±  -  XjUth  observation  is  tsicensored) ,  and 

(7) 

■  0  if  ZA  ■  Tj^ith  observation  is  censored) . 

We  assume  that  X^,  ...,  X^  are  independent  and  identically  distributed  (i.i.d.)  according 
to  the  continuous  failure  df  F,  Tj,  ...,  TR  are  i.i.d.  according  to  a  continuous  censoring 
df  H(x)  «  P(T  s  x),  and  the  T's  and  X's  are  mutually  independent.  The  censoring 
df  H  is  typically,  though  not  necessarily,  unknown  and  is  treated  as  a  nuisance  parameter. 

Model  (7)  is  useful  because  in  many  situations  the  data  are  analyzed  before  all  the 
subjects  have  experienced  the  end-point  event.  For  cocaagile,  in  situation  (i)  described 
earlier,  some  components  on  test  will  still  be  functioning  when  the  test  ends;  in 
situation  (ii)  some  patients  will  have  dropped  out  of  the  study  (by  moving  to  another 
city,  for  eranple)  and  others  will  still  be  alive  at  data-analysis  time;  in  situation 
(iii)  some  victims  will  still  be  awaiting  extrication  whan  the  experiment  ends. 
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In  Section  2  we  derive  our  test  statistic*  denoted  by  V°.  The  derivation  is 
similar  to  the  approach  of  HP (1975),  but  whereas  their  derivation  utilized  the  empirical 
df  to  estimate  F,  our  statistic  is  formed  by  using  the  Kaplan -Meier  estimator  (1958) 

£ KM  (1958) J  to  estimate  F. 

Using  results  of  Reid  (1979)*  we  show  (in  Appendix  I)  that  V0,  suitably  standardized* 
is  asymptotically  normal .  The  null  asymptotic  mean  of  V4"  is  0*  independent  of  the 

nuisance  parameters  v  and  H.  However*  the  null  asymptotic  variance  of  nV"  (given  by  (A.4)) 

*2 

does  depend  on  u  and  H  and  must  be  estimated  from  the  data.  A  consistent  estimator  tg  is 

2 

also  presented  in  Section  2.  (The  development  of  tq  is  given  in  Appendix  I.) 

In  Section  3*  by  comparing  the  statistic  with  the  KP(197S)  V*  statistic*  we  derive 
a  measure  of  the  efficiency  loss  incurred  due  to  the  presence  of  censoring. 

In  Section  4  we  apply  the  V*"  statistic  to  some  prostate  cancer  survival  data 
considered  earlier  by  Koziol  and  Green  (1976)  and  Hollander  and  Proschan  (1979) . 

2.  A  Test  for  Monotone  Mean  Residual  Life  Using  Randomly  Censored  Data 
Let  D(s*  t)  •  Rp(s)Rj.(t){cp(s)  -  Cp(t)}.  Then  D(s,  t)  is,  for  s  <  t,  a  weighted 
measure  of  the  deviation  from  Hq  (exponential  ity)  to  (drarl  alternatives) .  The 
weights  Rp(s)  and  Rp(t)  represent  the  proportions  of  the  population  still  alive  at  times 
s  and  t,  respectively.  Note  that  D(s,  t)  ■  0  for  all  s  s  t  if  and  only  if  Hq  is  true. 

HP (197 5)  were  thus  led  to  consider 

A(F)  -  /  /Rp(s)Rp(t){cF(s)  -  Cp(t)}dF(s)dF(t)  (8) 

as  an  average  value  of  the  deviation  D.  Straightforward  calculations  show  that  a(F)  can 
be  rewritten  as 

A(F)  -  |?{(%)F(u)  -  (3/2)F2(u)  ♦  (4/3)F3(u)  -  (l/3)F4(u))du.  (8*) 


m 


HP  (1975)  formed  their  statistic  by  replacing  F  by  the  eapirical  distribution  function  F 
in  (8). 

In  our  randomly  censored  model,  the  empirical  distribution  cannot  be  computed  but  we 

A 

can  compute  F,  the  ttl(19S8)  estimator  of  F.  Letting  -  0,  and 

Z(i)  <  Z(2)  <  ...  <  Z(n)  denote  the  ordered  Z's,  and  6^  the  6  corresponding  to  Z^,  we 

have 

L(x)  •  1  -  F(x)  -  n  C(n  -  i)(n  -  i  ♦  l)"1)  (i).  (9) 

h  (i:Z(i)sx> 

In  (9)  we  treat  Z^  as  an  uncensored  observation  whether  it  is  uncensored  or  not. 
Although  ties  have  probability  zero  under  our  assumptions,  in  practice  ties  will 
occur.  When  censored  observations  are  tied  with  uncensored  observations,  the  convention 
when  forming  the  list  of  ordered  Z^  's  is  to  treat  the  uncensored  members  of  the  tie  as 
preceding  the  censored  members  of  the  tie. 

KM (1958)  show  that  F  is  the  maximum  likelihood  estimator  of  F  in  the  censored 
nonparametric  model  where  no  parametric  assumptions  about  F  are  made.  Asymptotic 

A 

properties  of  F  are  studied  by  Efron  (1967),  Breslow  and  Crowley  (1974),  Meier  (1975), 
and  Peterson  (1977).  Peterson  also  presents  a  suraaary  of  properties  of  f . 

To  form  our  test  statistic,  we  replace  F  by  F  in  (8’) .  Since  A(F)  is  not  scale 
invariant,  in  order  to  make  our  test  scale  invariant  we  use  the  test  statistic 

V*  -  A(F)/C,  (10) 

m  m 

where  v  ■  ]Rp(x)dx  is,  under  the  assumption  that  the  mean  ^Rp(x)dx  is  finite  and  under 
suitable  regularity  on  the  amount  of  censoring,  a  consistent  estimator  of  u. 

For  computational  purposes,  use 

v  -  I  C  Jl\(n  -  ))/(n  -  )  *  l)j<0)3(Z(i)  -  *(l.1))-  dD 


A': 

:-*k 
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For  computational  purposes  it  is  also  helpful  to  rewrite  a(F)  as 

.  n  i-1  i-1  26m 

A(F)  -L  l  {-(1/6)  n  C(n  -  j)/(n  -  j  ♦  l)i  UJ  +  ft)  n  [(n  -  j)/(n  -  j  ♦  1)3  & 

i-1  j-1  j-1 

(12) 

i-1  46,... 

-0/3)  1K&1  -  j)/(n  -  j  ♦  l)j  •  (Z(i)  -  z(i.1)J. 

A  consistent  estimator  rl  of  the  null  asymptotic  variance  of  nV  is  developed 


in 


a2 


Appendix  I.  For  computational  purposes  tq  can  be  written  as 

t J  -  720"1  +  "j*n(n  -  i  +  D’V  -  i)‘1{72'1exp(-2i"1Z(i)) 


-  18*1exp  (-3i  *XZ  ^  )  +  M^eoqrt-Afi^Zpj) 

♦  45_iexp(-5i’1Z(i))  -  18‘1exp(-6i‘1Z(i)) 

♦  72"1exp(-8i*1Z(i))>  ♦  n{72'1exp(-2w'1Z(n))  (13) 

-  18‘1exp(-3M’1Z(n))  ♦  16’1«ocp(-4»i  ’XZ  ^ ) 

♦  45’1exp(-5G*1Z(n))  -  18”1e>cp(-6jl'1Z(n)) 

♦  72“1otp(-8il"1Z(n))}, 
where  u  is  defined  by  (11) . 

To  test  Hq  versus  at  the  approximate  o -level,  reject  Hq  in  favor  of  if 
nVtJ1  >  zo,  and  accept  Hq  otherwise.  Here  is  the  upper  a  percentile  point  of  a 
standard  normal  distribution. 

To  test  Hq  versus  HJ  at  the  approximate  a -level,  reject  Hq  in  favor  of  H£  if 
ifrvS-1  <  -ia,  and  accept  otherwise. 


Insisting  that  the  null  asymptotic  variance  of  nV\  as  given  by  (A.4)  of  Appendix  I, 
be  finite  imposes  a  restriction  on  the  amount  of  censoring  that  can  be  tolerated  by  the  test 
based  on  V°.  The  right -hand- side  of  (A.4)  fails  to  converge  if  R^Cx)  »  0{R^(x)>  as  x  ap¬ 
proaches  «.  Consider  this  condition  in  the  proportional  hazards  model  where  R^(x)  ■  (Pp(x)}Y. 
Here  y  is  viewed  as  a  censoring  parameter  since  the  probability  that  an  observation  will  be 
censored  is  P(«^  *  0)  ■  y/(y  +  1).  Then,  in  this  model,  the  right-hand-side  of  (A.4)  is 
finite  if  y  <  1.  Thus,  for  the  proportional  hazards  situation,  the  V0  test  will  be 
inappropriate  when  the  expected  proportion  of  censored  observations  is  greater  than  or 
equal  to  *j. 

3.  Efficiency  Lose  Due  to  Censoring 

Since  the  V°  statistic  introduced  in  this  paper  is  a  generalization  of  the 
HP (197 5)  V*  statistic,  we  find  it  natural  to  compare  the  power  of  the  V*  test  based  on  n 
observations  in  the  non-censored  case  with  the  power  of  the  V°  test  based  on  n*  observations 
in  the  randomly  censored  model. 

Thus,  let  F.  be  a  parametric  family  within  the  dmrl  class  with  F.  being  exoonential. 

9  0O 

Then  we  assume  the  randomly  censored  model  with  F  *  F  and  with  censoring  distribution  H. 

For  a  sequence  of  alternatives  ®n  *  +  c*1**  (with  c  >  0)  tending  to  the  null  hypothesis, 

let  0n(en)  be  the  power  of  the  approximate  a -level  V*  test  based  cm  n  observations  in  the 

uncensored  model,  and  let  Pn»(®n)  denote  the  power  of  the  approximate  a-level  V~  test 

based  on  n'  observations  in  the  randomly  censored  model.  Consider  n'  ■  h(n)  such  that 

lim  B  (ej  ■  lim  B  (6„) ,  where  the  limiting  value  is  between  0  and  1,  and  let 
n  n  n*  n 

k  *  lim  n/n*.  The  value  of  k  can  be  viewed  as  a  measure  of  the  efficiency  loss  due  to 
censoring.  The  value  of  k  is  adapted  from  Pitman's  (cf.  Noether,  1955)  measure  of  asymptotic 
relative  efficiency  but  the  interpretation  of  k  must  be  modified  because  the  tests  based  on 
V*  and  V0  arp  not  competing  tests  in  the  randomly  censored  model  (V*  cannot  be  applied 
to  the  data  arising  in  the  randomly  censored  model).  Roughly  speaking,  for  large 
n  and  dmrl  alternatives  close  to  the  null  hypothesis  of  exponentiality,  the  V°  test 
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requires  n/k  observations  from  the  randomly  censored  model  to  do  as  well  as  the  V*  test 
applied  to  n  observations  from  the  uncensored  model .  Since  V0  and  V*  have  the  same 
asymptotic  means,  it  can  be  shown  that  k  reduces  to  the  limiting  ratio  of  the  null 
asynptotic  variance  of  rS*  to  that  of  nV,  namely, 

k  de-f ‘  e^(V\  V*)  -  (210 'VtJ,  (14) 

2 

where  Tq  is  given  by  (A.4) .  Note  that  k  depends  only  on  the  censoring  distribution  H,  and 
not  on  the  parametric  family  Ffl  of  dmrl  alternatives.  Hence  we  use  the  notation  e^(Vc,  V*), 
rather  than  ep  H(VC,  V*),  in  (14). 

We  consider  the  cases  (i)  where  the  censoring  distribution  is  exponential,  (x)  *  1 
for  x  <  0,  R^(x)  *  exp(-Xx),  x  >  0,  and  (ii)  where  the  censoring  distribution  is 
piecewise  exponential,  R^(x)  *  1  for  x  <  0,  and  for  r  *  1,  ....  m,  Rj^(x)  *  c^expC-XjX) , 
sr_i  <  x  s  sr,  and  R^(x)  ■  c^pcpC-x^x)  for  sffl  <  x,  where 

cr  -  exp(-  l  Xi(sj.  -  si  X)  ♦  Xrsr  l)  and  sQ  »  0. 

X  1  o  s 

In  order  that  tq,  given  by  (A.4),  be  finite,  in  case  (i)  we  must  impose  the  restriction 
X  <  1.  Then  from  (A.4)  (with  w  -  1)  and  (14)  we  find 

e^O^,  V*)  -  (6/35) { (1  -  X)'1  -  6(2  -  X)"1  +  9(3  -  X)"1  ♦  4(4  -  X)'1 

(15) 

-  12(5  -  X)"1  +  4(7  -  X)'1)"1. 

Values  of  ^(V*”,  V*)  are  given  in  Table  1.  Note,  from  (15),  that  as  X  tends  to  0 
(corresponding  to  the  case  of  no  censoring),  as  expected  we  have  e^CV6,  V*)  tends  to  1. 

In  order  to  provide  a  reference  point  to  the  amount  of  censoring,  and  thereby 
facilitate  the  interpretation  of  e^CV0,  V*),  we  also  include  in  Table  1  the  value  of 
■  P(X  <  T)  ■  (1  +  x)  '1,  the  probability  of  obtaining  an  uncensored  observation  when  X 
is  exponential  with  scale  parameter  1  and  T  is  independent  of  X  and  has  the  censoring 
distribution  H^. 


When  the  censoring  distribution  is  again  in  order  that  Tq  be  finite  we  wust  inpose 
the  restriction  ^  <  1.  Straightforward  but  tedious  calculations  yield 

e^Or ,  V»)  -  (6/35) cd  -  y^expC-Cl  -  Xr)sr  l J-expC-(l  -  xr)sr3) 

-  6<2  -  -  VVl3'63^^2  •  xP*r3)  +  9<3  -  V"1^^3  •  Wl] 

-esxp[-(3  -  Xr)sr))  (16) 

♦  4(4  -  Xr)_1{expt-(4  -  XJ.)sr_1]  -exp[-(4  -  X^]}  -  12(5  -  XjJ'^expc- (5  -  Xr)sr_13 
-exp[-(5  -  xr)sr3)+4(7  -  Xr)'1{exp[-(7  -  xr)sr.1]-expC-(7  -  xr)sr)}]  _1, 


where  s, 


m+1 


Values  of  e^CV0,  V*)  are  also  given  in  Table  1.  Again,  as  a  reference  point  for 
the  amount  of  censoring  under  we  include  in  Table  1  values  of  ■  P(X  <  T)  when  X 
is  exponential  with  scale  parameter  1,  and  T  is  independent  of  X  with  distribution 
Direct  calculations  show 

hh-1  , 

•  1  -  lcrx Cexp(-(xr  +  l)sr  l)-exp{-(xr  ♦  l)sr>]. 


where  s. 


m+1 


(Insert  Table  1  here.J 


4.  Exempt* 

The  data  in  Table  2  were  analyzed  in  Hollander  and  Proschan  (1979)  [HP(1979)J  and 
are  an  updated  version  of  data  considered  by  Koziol  and  Green  (1976) .  The  data  corre¬ 
spond  to  211  state  IV  prostate  cancer  patients  treated  with  estrogen  in  a  Veterans 
Administration  Cooperative  Urological  Research  Grow  (1967)  study.  At  the  March,  1977 


it.  <f  % 


closing  date  there  were  90  patients  who  died  of  prostate  cancer,  105  who  died  of  other 
diseases,  and  16  still  alive.  The  latter  121  observations  are  treated  as  censored  obser¬ 
vations  (withdrawals). 

Figure  1  of  HP(1979)  gives  a  plot  of  the  KM  survival  function  Rp  for  the  data  of 
Table  2.  H>(1979)  also  applied  several  goodness -of- fit  tests,  including  ones  proposed 
by  Kbziol  and  Green  (1976),  Hyde  (1977)  and  a  test  developed  in  H>(1979),  to  test  the 
null  hypothesis  that  F  is  exponential  with  mean  100  months.  (As  stated  by  Kbziol  and 
Green  (1976),  prior  experience  suggested  that  had  the  patients  not  been  treated  with 
estrogen,  their  survival  distribution  for  deaths  from  cancer  of  the  prostate  could  be 
taken  to  be  exponential  with  mean  100.) 

Although  the  tests  applied  in  HP(1979)  tended  to  support  the  postulated  exponential 
distribution  with  mean  100,  the  two-sided  P  value  for  the  Koziol -Green  test  was  .14 
suggesting  that  perhaps  a  different  model  could  be  more  appropriate.  Furthermore,  Gregory 
(1979)  analyzes  the  Table  2  data  by  new  goodness -of-fit  tests  (which  will,  for  some  alter¬ 
natives,  be  more  powerful  than  those  proposed  by  JCoziol  and  Green  (1976),  Hyde  (1977),  and 
HP  (1979))  and  his  tests  do  not  support  the  postulated  exponential  distribution  with  mean 
100. 

Figure  1  of  this  paper  plots  the  Kaplan -Meier  analogue  ?f(0  of  the  empirical  mrl 
function  (defined  by  replacing  Sp  by  Rp  in  (2))  for  the  data  of  Table  2.  Specifically 

*F(x)  •  C(f  Rp(u)du)/Rp(x))j«ICx<Z^  j. 

Note  that  Ip,  for  the  data  of  Table  2,  tends  to  decrease  up  to  around  25  months,  then  tends 
to  increase  up  to  about  70  months,  and  then  decreases  again. 

Figure  1  suggests  "wearout",  but  does  not  strongly  suggest  a  etrietly  decreasing 
mean  residual  life. 

Application  of  the  V6  test  tends  to  confirm  the  "eye-ball”  impression  one  gets  from 
Figure  1.  We  obtain,  from  (10)  -  (12),  Vs  •  .027,  and  from  (13),  ■  .066.  Ns  then  find 
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(211)1VcTq1  •  1.S2  with  a  corresponding  one-sided  P  value  of  .064.  Thus,  with  this  ob¬ 
jective  analysis,  we  find  that  the  test  suggests  "wearout"  in  the  tarl  direction,  but  the 
test  does  not  strongly  suggest  a  strictly  decreasing  mean  residual  life. 

[Insert  Table  2  here.] 

[Insert  Figure  1  here.] 
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Appendix  I 

A  Consistent  Estimator  of  the  Null  Asymptotic  Variance  of  V6 
Using  straightforward  algebra  and  integration  we  can  write 

^  -(ACF)!*"1)  -  GO"1  JxJ(F(x))dF(x),  CA.1) 

where 

J(u)  -  -a/2)  -  ACFJy"1  ♦  3u-  4u2  ♦  (4/3)tt3. 

The  tern  JxJ(F(x))dF(x)  is  a  Kaplan-Meier  "L-statistic"  and  results  of  Reid  (1979)  can  be 
used  to  show 

n*  JxJ(F(x))dF(x)  -*•  N(0,o2)  CA.2) 

where 

-  XAV  . 

oZ  -  //RpCxjRpOrpCFWjjCFCy))  ^  {R^CuJRpMJ^dFCuJdxdy,  (A.3) 

where  RjOO  •  RpCuJRjjCu)  and  x  a  y  -  nin(x,y) .  Thus  from  (A.1)  -  (A.3)  and  Slutsky's  S  * 
theorem  (cf.  Cramdr,  1946,  pp.  254-255),  we  have 

n^CV*"  -  (ACPjw"1})  N(0,o"2o2). 

Under  1^,  it  can  be  shown  that  y"Zo2  reduces  to  the  right-hand-side  of  (A.4),  namely 

***•  <M) 

where 

g(z)  -  ze^tz  -  6*2  ♦  9z3  ♦  4z4  -  12z5  ♦  4z7). 

When  there  is  no  censoring,  i.e.  R^(x)  s  1,  0  <  x  <  •,  the  right-hand-side  of  (A.4)  re&ices 
to  1/210,  agreeing  with  the  result  for  the  asymptotic  variance  of  nV*,  the  ff&975) 
statistic  designed  for  the  uncensored  case. 

A  consistent  estimator  t q  o f  Tq  is  obtained  by  replacing  Rg,  in  (A.4),  by  Rg»  the 

empirical  survival  function  of  the  Z*s;  Rg(x)  -  n’1  •  {number  of  Z's  >  x).  The  resulting 

estimator  is  given  by  the  right-hand-side  of  (W)» 
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TABLE  1 

Efficiency  Losses  Under  Exponential  (H^)  and  Piecewise 
Exponential  (Hj)  Censoring 


«2 

x: 

1/2 

1/3 

1/4 

1/10 

.248 

.435 

.552 

.805 

V 

.667 

.750 

.800 

.909 

«2 

n  •  1  Sj: 

1 

1 

1 

2 

(1/2,1/S) 

(1/3.1/2) 

(1/3, 1/4) 

(1/2, 1/3) 

e^flf.V*): 

.374 

.290 

.513 

.326 

V 

.685 

.728 

.763 

.671 

m  -  2  (SpS2): 

(1/2,1) 

d/2,1) 

(1/2,1) 

(*1»*2**3^  * 

(1/2, 1/3, 1/4) 

(1/4, 1/3, 1/2) 

(1/2,1/5,1/10) 

e^C^.V*): 

.475 

.302 

.645 

■V 

.718 

.750 

.765 

fh'. 


TABLE  2 


Survival  Tines  and  Withdrawal  Tines  in  Months  for  211  Patients 
(with  number  of  ties  given  in  parentheses) 


&nrvival  tines:  0(3),  2,  3,  4,  6,  7(2),  8,  9(2),  11(3),  12(3),  15(2),  16(3),  17(2),  18, 
19(2),  20,  21,  22(2),  23,  24,  25(2),  26(3),  27(2),  28(2),  29(2),  30,  31,  32(3),  33(2),  34, 
35,  36,  37(2),  38,  40,  41(2),  42(2),  43,  45(3),  46,  47(2),  48(2),  51,  53(2),  54(2),  57,  60, 
61,  62(2),  67,  69,  87,  97(2),  100,  145,  158. 

Withdrawal  tines:  0(6),  1(5),  2(4),  3(3),  4,  6(5),  7(5),  8,  9(2),  10,  11,  12(3),  13(3), 
14(2),  15(2),  16,  17(2),  18(2),  19(3),  21,  23,  25,  27,  28,  31,  32,  34,  35,  37,  38(4),  39(2), 
44(3),  46,  47,  48,  49,  50,  53(2),  SS,  56,  59,  61,  62,  65,  66(2),  72(2),  74,  78,  79,  81,  89, 
93,  99,  102,  104(2),  106,  109,  119(2),  125,  127,  129,  131,  133(2),  135,  136(2),  138,  141, 
142,  143,  144,  148,  160,  164(3). 
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